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An improvement of pointing precision of the DSS-14 antenna is planned for the 
near future. In order to analyze the improvement limits and to design new con- 
trollers, a precise model of the antenna and the servo is developed, including a finite 
element model of the antenna structure and detailed models of the hydraulic drives 
and electronic parts. The DSS-14 antenna control system has two modes of opera- 
tion: computer mode and precision mode. The principal goal of this investigation 
is to develop the model of the computer mode and to evaluate its performance. 
The DSS-14 antenna computer model consists of the antenna structure and drives 
in azimuth and elevation. For this model, the position servo loop is derived, and 
simulations of the closed-loop antenna dynamics are presented. The model is sig- 
nificantly different from that for the 34-m beam-waveguide antennas. 


I. Introduction 

The DSS-14 antenna control system model consists of the antenna structure, antenna drives in azimuth 
and elevation, and the position servo loop. Each drive, in turn, consists of gearboxes, hydraulic servo 
(active and passive valves, hydraulic lines, and hydraulic motors), and electronics boards (amplifiers 
and filters). The DSS-14 antenna control system model was developed by R. E. Hill [1,2]. In the present 
development, we obtain a more precise model that allows for accurate simulations of the antenna pointing 
errors and allows simulation of the intermediate variables, such as torques, currents, wheel rates, truss 
stresses, etc. We incorporate the finite element structural model with free rotation in azimuth and 
elevation, in a manner similar to the 34-m antenna models [3-5], that involves cross-coupling effects 
between azimuth and elevation, wind pressure on the dish, and pointing error model. The hydraulic part 
involves a recent development in modeling of the hydraulic components by R. Bartos [6-8]. 

The rate loop model consists of the elevation and azimuth drives and the antenna structure. Each drive 
consists of three major components: the electronics boards, hydraulic system, and gearbox. A model of 
each component is derived separately, then put together, forming the drive and rate loop models. Finally, 
the position loop is closed to obtain the position loop model. 


II. Drive Model 

A block diagram of the drive model is shown in Fig. 1, where N t is the ratio between motor rate 
and tachometer (pinion) rate; r, rad/s, is the rate input to the drive; i, A, is the hydraulic active valve 
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Fig. 1 . Block diagram of the antenna drive. 


solenoid current; T 0 and T, N-m or lb in., are the gearbox and on-axis torques, respectively; and d tac h and 
6 m , rad/s, are the pinion and motor rates, respectively. The state-space representation of the electronic 
board, hydraulic system, and gearbox are derived in the following sections. 

A. Electronic Board 

A schematic diagram for the electronic board is shown in Fig. 2. The inputs are the rate command 
r, rad/s, and the tachometer rate 6 tac h , rad/s. The output is the solenoid valve current i. The scaling 
factors, k r and k t , convert the inputs into the command voltage, v r , and tachometer voltage, v t . The 
subsystem, G t , is the tachometer circuit: it transforms the tachometer voltage, vt, into the voltage, v t0 . 
The subsystems with the transfer functions G r i and G r 2 are the rate amplifier circuits: they transform 
the command voltage, v r , and the tachometer voltage, v t0 , into the error voltage, v s . The subsystem with 
the transfer function, G s , is the valve driver amplifier circuit, with the error voltage, v 3 , as the input and 
the valve current, i, as the output. 



Fig. 2. Block diagram of the electronic board. 


The following transfer functions of each of the four components are derived in the Appendix. The 
transfer function Gto for azimuth is 


G to = 0.151 


( 1 ) 


and, for elevation, is 


G to = 0.127 


( 2 ) 


The transfer functions G rl (from v r to v s ) and G r2 (from v to to v 3 ) are 
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G r 1 = 6.20 Go ' 
G r2 = — 4.65 G 0 t 


where 


1 4- 0.400s 
1 4- 4.205s 


is the transfer function of a lag compensator. The transfer function G a is 


G s = 4.42 x 10~ 5 


( 3 ) 


( 4 ) 


( 5 ) 


The scaling factors, k r and k t , are k r = 1212.6 V/rad/s and kt = 2.5 V/rad/s. Thus, the command 
transfer function from the rate command r to the solenoid current i s is 

G r = k r G rl G s = 0.3323 G 0 (6) 

where G 0 is defined in Eq. (4). The tachometer transfer function from the tachometer rate Otach to the 
solenoid current i s is 


G t = k t G to G r2 G s (7) 

q _ ( —0.7750 x 10~ 4 G o for azimuth /g.. 

4 \ —0.6525 x 10~ 4 G o for elevation 

In order to check the correctness of the derivation, note that the ratio G r /G t should be equal to N 0 , where 
N 0 is the tachometer-to-axis ratio ( N 0 = 4287.5 for azimuth and N a = 5083.6 for elevation). Indeed, 
from Eqs. (6) and (7), one obtains 


G r ( —4287.7 = N 0 for azimuth 

G t \ —5092.7 = N 0 for elevation 


Finally, the state-space representations of the transfer functions G r and Gt (for azimuth and elevation) 
are easily obtained with the standard Matlab command in the form 


Xfo — AfoXb 4 “ B V 4 “ Bb2@tach 

> 

i = GfoX 4 “ T Db 2 6tach t 


( 10 ) 


The plot of the transfer function in azimuth (magnitude and phase) from r to i is shown in Fig. 3. The 
transfer function for elevation is identical. The plots of the transfer functions in azimuth and elevation 
from Otach to i are shown in Fig. 4. The magnitudes drop in the frequency range from 0.01 to 0.1 Hz due 
to implementation of the filter G c . 
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Fig. 3. Magnitude of the electronic board transfer function from the rate input to the 

board current. 



B. Hydraulic System 

The hydraulic servo system model was presented by R. E. Hill in [1] and [2]. Here we take a different 
approach, based on the recent investigations of hydraulic components by R. Bartos (see [6-8]). A block 
diagram of the DSS-14 hydraulic system is shown in Fig. 5. It consists of the hydraulic motor, shorting 
valve, hydraulic lines A and B, passive servo valves, and active servo valves. It has two inputs, servo 
valve current i and motor rate 6 m , and one output, motor torque T a . The equations for each component 
are derived separately based on the work of Bartos [6-8]. Basically, these models are nonlinear ones; 
however, we linearize them in order to model the antenna linear regime of operation. 

1. Active Servo Valve. This valve model has the input, i, A, and two outputs, q av — the flow rate 
out of port a, cm 1 * 3 /s, or in. 3 /s, and qb v , the flow rate out of port b, cm 3 /s, or in. 3 /s. From [6], one obtains 

Qav ‘^Co^oQav T U 0 ([ av — Ld 0 k a i (11a) 

Qbv “ Qav (lib) 


where ( 0 = 0.8 is the damping ratio, ui 0 = 345.6 rad/s is the valve natural frequency, and k a — 59,200 
— 97,300 cm 3 /s/A (23,300-38300 in. 3 /s/A) is the valve gain. The lower value is the Bartos estimate, 
while the upper value is the Hill estimate [2]. The values of the parameters are listed in Table 1. 
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Fig. 5. Block diagram of the hydraulic drive. 


Table 1. Parameters of the active servo valve (line A and line B). 


Drive 

Co 

O>o, 

rad/s 

fta.5 

cm 3 /s/A 

in. 3 /s/A 

Azimuth 

0.6 to 0.8 

345.56 

59,200 to 97,300 

23,300 to 38,300 

Elevation 

0.6 to 0.8 

345.56 

59,200 to 97,300 

23,300 to 38,300 


Introducing the new variable q v = q av — qb v , one obtains from Eq. (11) 

q v + 2( 0 u 0 q v -I ~ = 2w 0 fc 0 t (12) 

The differential variable q v and the other differential variables introduced allow one to further simplify 
the analysis without loss of accuracy and to get rid of the “parasitic” variables, such as tank pressure, 
supply pressure, and case pressure. 

2. Shorting Valve. The pressures p a and pb , kPa (lb/in. 2 ), are the inputs to the shorting valve, and 
the flows q as and qbs, cm 3 /s (in. 3 /s), are its output (see Fig. 5). The linearized relationship between the 
inputs and outputs is as follows: 


Qas — ks(Pa Pb) 1 


Qbs — Qas 


(13) 


where k s is the valve gain, k s = 0.0007 — 0.007 cm 3 /s/kPa (0.0003-0.003 in. 3 /s/psi), both in azimuth and 
elevation. 
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Introducing new differential variables p = p a ~Pb and q s = q as — qbs, one obtains Eq. (13) in the form 


q s = 2k s p 


(14) 


3. Passive Servo Valve. This valve has four inputs: pressures p a and pb, supply pressure p s , and 
tank pressure p t . The last two are supplementary constant inputs that can be removed from the analysis. 
The valve has two outputs: flows q ap and qb p . Its linearized input-output relationship is as follows: 


Qap = k p i(p a - p s ) + k p2 {p a - Pt) (15a) 

qbp = k p i(p b - p s ) + k p2 (pb - Pt) (15b) 


where the gains are k p \ = 0.0055 cm 3 /kPa (0.00233 in. 3 /psi) and k p2 = k p i, the supply pressure is 
17,240 kPa (2500 psi), and the tank pressure is 345 kPa (50 psi). These values are identical for azimuth 
and elevation. 

Introducing q p = q ap — qbp , and recalling that p = p a ~Pb , one obtains Eqs. (15a) and (15b) as follows: 


q P = (k p i + k p2 )p — 2 k p p 


(16) 


where, for simplicity of notation, we denote k p — k p i — k p 2 ■ 

4. Hydraulic Motor. The motor is described in [8]. From Fig. 5, it follows that the motor has four 
inputs and four outputs. The inputs are pressures p a and pb, case pressure p c , and motor rate 6 m , rad/s. 
The outputs are flows q a and qb, leakage to the case q c , and motor torque T 0 , N-m (or lb in.). Following 
[8], one obtains the flow q a from Eq. (59) of [8]: 


q a = q a i + q a 2 + q a 3 


(17) 


But, from Eq. (40) of [8], 


Qal — DQ m 


(18) 


where D = 6.3 cm 3 /rad (0.3836 in. 3 /rad) is the motor “displacement.” From Eq. (52) of [8], one obtains 


?a2 = ~~~~{Pa ~ Pc) (19) 

A 4 

where k a2 = from 6.35 x 10 -5 to 18.3 x 10 -5 cm 3 (from 2.5 x 10 -5 to 7.2 x 10 -5 in. 3 ) is the leakage 
constant (assumed to be 10 -4 cm 3 , or 4 x 10 -5 in. 3 ); p — from 2.8 x 10 -4 to 2.8 x 10~ 3 kPa s (from 
4 x 10 -5 to 4 x 10 -6 lb s/in. 2 ) is the absolute fluid viscosity (assumed to be 10 -4 cm 3 , or 4 x 10 -5 in. 3 ); 
and p c = 221 kPa (32 psi). From Eq. (58) of [8], one obtains the leakage from port A to port B, q a 3: 


ka3 ( , 

qa 3 = \Pa Pb) 

A* 


(20) 
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where k a 3 is the constant of proportionality determined through experiments. It is assumed to be equal 
to fc a 2 ) k/i '3 = k a 2 * 


Combining Eqs. (17) through (20), one obtains 


i k a 2 T kaZ k a 3 k a 2 

q a = D6 m + p a p b p c 

p p p 


( 21 ) 


From [8], Eq. (60), one obtains the flow rate <?{,: 


96 = 9bl + 962 “ 9 o3 


( 22 ) 


It follows from [8], Eq. (41), that 


96i — 9ai — D9 n 


(23) 


and from [8], Eq. (53), that 


962 = —(Pb-Pc) 
t 1 


(24) 


Combining Eqs. (22), (23), (24), and (20), one obtains 


A k a 3 k a 3 'r kb2 h,2 

q b = -DO p a + pb p c 

p p p 


(25) 


The motor torque T a is obtained from Eq. (28) of [8] by neglecting Coulomb friction and inertia torques 
(the latter are included in the gearbox model): 


T 0 = T P + T f (26) 

where T p is the torque generated by the motor and Tj is the viscous friction torque. The linearized 
Eq. (10) of [8] gives the torque generated by the motor: 

Tp — D(Pa ~ Pb) ( 27 ) 

and from Eq. (24) of [8], one obtains the viscous friction torque: 

T f = —k v Dp9 (28) 


where k v — 0.0438 is a dimensionless viscous friction coefficient. Combining Eqs. (26) through (28), one 
obtains 


T 0 = Dp a - Dp b - k v Dp9 


(29) 
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Define q = q a — qb\ then from Eqs. (21) and (25), one obtains 


o no i *^<*2 

q = 2 D9 m + p 

M 

T 0 = Dp - k v Dp9 m 


( 30 ) 


The motor parameters are given in Table 2. 


Table 2. Parameters of the hydraulic motor (line A and line B). 


Drive 

fi, kPa s 

H, lb s/in. 2 

D, 

cm 3 /rad 

D, 

in. 3 /rad 

k-y 

fc Q 2 , cm 3 

k a 2 , in. 3 

Pc, 

kPa 

Pci 

psi 

Azimuth 

2.8 X 10~ 4 

0.4 x 10~ 6 

25.2 

1.52 

0.0438 

4.1 x 10“ 4 

2.5 x 10- 5 

220 

32 


to 

to 




to 

to 




28 x 10" 4 

4 x 10" 6 




11.8 x 10“ 4 

7.2 x 10“ 5 



Elevation 

2.8 x 10~ 4 

0.4 x 10~ 6 

25.2 

1.52 

0.0438 

4.1 x 10“ 4 

2.5 x 10- 5 

220 

32 


to 

to 




to 

to 




28 x 10- 4 

4 x 10~ 6 




11.8 x 10- 4 

7.2 x 10- 5 




5. Hydraulic Line. There are two lines: A and B. A model for line A is developed, and the model 
for line B is similar (index “a” should be replaced with ,s b”). Line A has four inputs, flows q a , q a v,q a pv, 
and q asv , and a single output, pressure p a (refer to Fig. 5). From [7], one obtains the line- A model as an 
integrator, with the negative feedback signs as in Fig. 5: 


Pa — k[ a ( q a -f- (Jav qap i/as) 


(31a) 


Similarly, the line-B model is obtained: 

Pb = kib{ qb qbv — qbp qbs} 


(31b) 


As before, defining p = p a — Pb, one obtains 


p = ki{-q + q v - q p - q s ) 


(32) 


In these equations, the gains are 



(33) 


where /3 is the effective bulk modulus (capacitance of the line), /? = 1.29 x 10 6 kPa (1.87 x 10 5 psi), and 
v Q is the total volume, v Q — 27, 200 cm 3 (1660 in. 3 ), so that ki — 47.3 kPa/cm 3 (113 psi/in. 3 ). The values 
are collected in Table 3. 
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Table 3. Parameters of line A and line B. 


Drive 

/3, kPa 

13, psi 

v 0 , cm 3 

v 0t in. 3 

Line A 

Azimuth 

Elevation 

1.29 x 10 6 
1.29 x 10 6 

1.87 x 10 5 
1.87 x 10 5 

27,200 

24,500 

1653.6 

1498.1 

Line B 

Azimuth 

Elevation 

1.29 x 10 6 
1.29 x 10 6 

1.87 x 10 5 
1.87 x 10 5 

27,700 

24,600 

1690.4 

1499.9 


6. Hydraulic System Model. The model of the hydraulic system is derived by combining its 
elements (active servo valve, shorting valve, passive valve, hydraulic lines, and hydraulic motors). By 
introducing the new differential variables, the block diagram in Fig. 5 is simplified to the one in Fig. 6. 
A detailed block diagram of the hydraulic system is shown in Fig. 7. Combining Eqs. (12), (14), (16), 
(32), and (30) (or, alternatively, using the block diagram in Fig. 7), and defining the new state vector 
x h = [xi,X 2 ,X 3 ] t , with three states, xj = q v ,x 2 = q v , and X3 = p, and defining the input current, i, 
motor rate 6 m , and the single-output motor torque, T 0 , one obtains 


Xh — A/,.7'/, + Bho@m T ^ hi i 
7' 0 = Ch^h Dho@m T Dhi j 


(34a) 


where 



Fig. 6. Simplified block diagram of the hydraulic drive. 
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Fig. 7. Detailed block diagram of the hydraulic system. 



(34b) 


The plots of the magnitudes of the transfer function in azimuth and elevation from i to T a are shown 
in Fig. 8. The plots of the magnitudes of the transfer functions in azimuth and elevation from 9 m to T a 
are shown in Fig. 9. 

C. Gearbox Mode! 

The gearbox model was described in detail in [5], and its block diagram is given in Fig. 10. In this 
diagram, T 0 is the motor torque, 9 P is the antenna angular rate, w m is the motor rate, T is the gearbox 
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Fig. 8. Magnitude of the hydraulic drive transfer function from the solenoid current to 

the output torque. 



Fig. 9. Magnitude of the hydraulic drive transfer function from the motor rate to the 

output torque. 



Fig. 10. Block diagram of the gearbox. 


torque, J m is the motor inertia, k g is the gearbox (output) stiffness, and N is the gearbox ratio. This 
model has two inputs, the motor torque, T c , and the wheel (pinion) angular rate, 8 P , and a single output, 
the gearbox torque, T. 

The equations for this system are as follows: 
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Jm&m — T 0 


N 


(35a) 


Denoting the state variables xi = u> m and x 2 = T, one obtains 

• - ~ X2 i ll. 
Xl " NJ m + J m 


(35b) 


(36a) 


± 2 


kgX 1 

“j\T 


kgOp 


(36b) 


Defining the gearbox state as x g = [xi X 2 ] T , input T e and 0 P , and output T and uj m , 
gearbox state-space representation (A g . B g ,C g ): 

Xg = AgXg + B gl T 0 + B g 2 0 p 

T = CglXg > 

@m = Cg2 x g . 


one obtains the 


(37a) 


where 


Ag 



-1 -| 
NJm 
0 




> 


< 7 * 1=10 1 ] 
< 7*2 =[1 0 ] 


(37b) 


D. Drive Model 

The drive model is obtained by combining the state-space representation of the electronic board, 
Eq. (10); the hydraulic system, Eq. (34); and the gearbox, Eq. (37), according to the block diagram in 
Fig. 1. Defining the drive state vector Xd = [xj ,x^,xJ] T , we obtain the state equations 
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> 


(38a) 


id = A d x,i + Bd r r + BdtOp 

T = CdXd 


where 


A d = 


A b 

BhiCb 


0 

Ah 


BiaCgi 


BhoC'g2 F 


N t 

BhiDbeCgi 


N t 


B d r 


B d t 


B g iD h iCb B g iC h A g + B gl D ho Cg 2 + BalD ^ b2C ^ 


Bbi 

BhiDbi 

L B gl D hi D bl 


0 

0 

L B g2 


C d = [0 0 C gl ] 


(38b) 


The plots of the magnitudes of the transfer function in azimuth and elevation from r to T are shown 
in Fig. 11. The plots of the transfer functions in azimuth and elevation from 9 P to T are shown in Fig. 12. 


III. Structure Model 

The structural model is derived from the finite element model of the antenna structure with free 
rotations with respect to the elevation and azimuth axes. The finite element model consists of the 
diagonal modal mass M m (p x p), diagonal natural frequencies matrix f l(p x p), diagonal modal damping 
matrix Z(p x p), and modal matrix $(m x p),p < m, which consists of p eigenvectors fa (mode shapes), 
i = 1, 




(39) 


Let the finite element model have m degrees of freedom, with s inputs u(t), where uissxl vector, 
and with r outputs y(t), where y is r x 1 vector. If the input matrix is B 0 (m x s), the output matrix 
for displacement is C oq (r x m), and the output matrix for rates is C ov (r x m), then the input-output 
relationship is given by the following second-order differential equation: 


Qm T 22fh) m F Qrn — ltl $ B 0 u 

> 

Vm ~^oq^9m F C ov ^Q m 


(40) 


Define the state variable x as follows: 
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Fig. 1 1 . Magnitude of the drive transfer function from the rate command to 
the output torque. 



Xi 


Qm 

_ x 2 . 


Qm 


(41) 


where q m and q m are modal displacements and rates (such that q = $g m ;g is the actual displacement); 
then Eq. (40) can be presented as a set of first-order equations: 


±i = x 2 


x 2 = - fl 2 xi - 2Z£lx 2 + M^QBoUs > 
Vs i C ov (j)X2 ^ 


or in the following form: 


X s AgXg BgUg 

> 

Vs —C s x 


(42) 


(43a) 
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where 


A s = 


0 

-ft 2 


I 

-2 Zft 


B. 


M-WBo 


C s =[C oq $ C ov $] 


(43b) 


is the sought state-space model in modal coordinates. In our case, u s = [T a T e ], where T a and T e 
are torques at azimuth wheels and elevation pinions, respectively. The structure output consists of the 
elevation and azimuth encoder angles and rates, pinion angles, elevation and cross-elevation pointing 
errors, and other structural variables of interest. Two outputs, 9 pa and 9 pe , the pinion rates in azimuth 
and elevation, are of special interest. Thus, the structural state-space equations are as follows: 


x s — A s x s + B sa T a + B se T e 
9pa = (-'pa X s 

8 p e — ('p a X s 

y = C s x j 


(44) 


The modal data obtained from the finite element model consist of 150 natural frequencies, w*; modes, 
<j>i\ and modal masses, m m j, i = 1, • • • , 150. Additionally, based on the measurements, the modal damping 
is assumed to be 1 percent, i.e., Q = 0.01. Based on this information, the state matrix A s , as in Eq. (43b), 
is determined by introducing the matrix of natural frequencies, ft = diag(uJi), and modal damping, 
Z = diag(£i), i — 1, • • ■ , 150. 

The determination of matrices B s and C s is presented here for the azimuth wheel torque input and 
the azimuth wheel rate output. For the azimuth wheel torque input, consider the azimuth wheel of radius 
r a and the azimuth rail of radius R a - Let nodes n i be located at the contact point of the wheel and the 
rail. The torque applied to the wheel generates the force F a at node n\. The force is tangential to the 
azimuth rail. Assuming a rigid pinion, the force F a applied to the wheel is 



(45) 


This force has x and y components, F ax and Fay [see Fig. 13(a)], such that 


X a 

Fax F a COS Oi a — — COS QIq 

r a 


T 

Fay = F a sin a a = — sin a a 

r a 


> 


(46) 
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Fig. 13. Forces and rates at the azimuth pinion: (a) forces and (b) rates. 

and a a is the angle marked in this figure. Let e x and e y denote the unit vector (all but one component 
are zero, and the nonzero component is equal to one), with the unit component at the location of the x 
and y displacement of node n\ in the finite element model. The input, F, to the finite element model 
is F — F ax e x + F ay e y . Therefore, B 0 follows from the decomposition of F, such that F = B a T a . Prom 
Eq. (46), it follows that 


B 0 = cos a a + — sin a a 

r a r„ 


(47) 


Next, from Eq. (43b), it follows that the nonzero (lower) part of B, after introduction of Eq. (47), is 


M l <b T B a = - M 1 cosa a + M* ' '-' y - sin a„ 


_ 1 $ T e x j $ T e y . 

one I A/f ^ oi' 


M m l — cos a a + — sin a 0 

T a T a 


(48a) 

(48b) 


where <j) x and <j> y are vectors of modal components of x and y displacements at node n i : 


0x — ^ Cx — < t > x2i ' ' ' i 0xl5o] 


T \ 


> 


(49) 


^ e y ~ [tfiyh 4 > y2) * * * > 0yl5o] J 

where <\> xt and 4> yi are x and y displacements of mode i at node ni. Therefore, from Eqs. (43b) and (48b), 
one obtains 


B s 


0 

~M m l — cos a a + M" 1 ^ sin a a 
r a r a 


(50) 


The output matrix derivation is presented here for the wheel rate, 6 pa . The wheel rotation is 
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(51) 



Va 

r a 


where v a is the tangential velocity of the wheel at the contact point [see Fig. 13(b)]. If v x and v y are x 
and y components of v a , and a a is the angle marked in this figure, then 


v a = —v x cos a a + Vy sin a a 


therefore, 



— — cos a a + 
r n . 



Q 


(52) 


(53a) 


and in modal coordinates 


^pa 


$el L . V 

— cos ot a H sin a a q m = 

r a r a / 


4>t <t>y \ 

cosa a + — Sin a a q m 

r a r a / 


(53b) 


Finally, the matrix C s , according to Eqs. (43b) and (53b), is 


C s 


4>l 

0 — — cos a a 

TV 


sma„ 


(54) 


The structural model consists of m = 150 modes or 300 states. Modes not participating in system 
dynamics are eliminated. Observability and controllability properties in the balanced representation are 
used to determine insi g nifi cant, modes. The balanced representation [9] is a state-space representation 
with equally controllable and observable states. The Hankel singular value is a measure of the joint 
controllability and observability of each balanced state variable. The stated with small Hankel singular 
values are deleted as weakly excited and weakly observed, causing minimal modeling error. 

For flexible structures with small damping and distinct poles, the modal representation is almost 
balanced, c.f. [10-12], and each mode is considered for the reduction separately. For a structure with 
m modes, matrix B s has 2 m rows, and C s has 2m columns. Denote b s as the last m rows of B s , c q as 
the first m columns of C s , and cv as the last m columns of C s . Then b si is the ith row of b s ,c q i is the 
zth column of c q , and tv* is the ith column of cv. Denote /3^ = & si 6T,a: gi = <? qi c q *, and a ri — c^c ri . The 
Hankel singular value for the ith mode is given in [11] and [12]: 


it 


WbiPsi y/rfiP%i + w li^ 2 i a li 


40? 


(55) 


where the weighting factors wu > 0, w qi > 0, w r i > 0, and i = 1, ■ • • , m. 

Care should be taken when determining Hankel singular values. Units should be consistent; other- 
wise, some inputs or outputs receive more weight in Hankel singular-value determination than necessary. 
Consider, for example, the azimuth encoder reading in arcseconds and the elevation encoder reading in 
degrees. For the same angle, the numerical reading of the azimuth encoder is 3600 larger than the eleva- 
tion encoder reading; hence, the elements for the azimuth output are much larger than those for elevation. 


129 



On the other hand, some variables need more attention than others: Pointing error and encoder readings 
are the most important factors in the antenna performance; hence, their importance has to be emphasized 
in mode evaluation. For consistency of units and importance of variables, the weighting factors Wbi,w qi , 
and w r i are introduced. Typically, weights are set to 1. 

For each mode, the Hankel singular value is determined and used to decide on the number of modes in 
the reduced structural model. For the rigid body modes, Hankel singular values tend to infini ty; hence, 
rigid body modes are always included in the reduced model. Hankel singular values of the 150 modes of 
the antenna model are plotted in Fig. 14. The reduced order model consists of 24 modes: 2 rigid-body 
modes and 22 flexible modes. 

The plots of the transfer function in azimuth and elevation (magnitude and phase) from the wheel 
(pinion) torque T to the axis rate 6 are shown in Fig. 15. They show that the azimuth transfer function 
has low frequency resonances (about 1.2 and 2.2 Hz), which are absent in the elevation transfer function. 




Fig. 15. Magnitude of the transfer function of the antenna structure: (a) direct coupling 

and (b) cross-coupling. 
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SV, Rate Loop Model 

A rate loop block diagram is presented in Fig. 16, where T a and T e denote the drive torques, 0 pa and 
6 pe denote pinion rates, and r a and r e are rate commands in azimuth and elevation, respectively. The 
state-space equations are combined from the state equations of the azimuth and elevation drives [see 
Eq. (38) and add subscript “a” for the azimuth drive and subscript “e” for the elevation drive] and the 
structure [see Eq. (44)]. Combining them, and defining the rate loop state vector x r as x r = [xd a , x de ,x s }, 
where Xda and Xd e are azimuth and elevation drive states, one obtains the rate-loop state-space equations: 


Xr = A r X r + B ra'T a + B 

re T e 


y — C r x r 

► 

0 a = CdXf 

Q e = C e x r j 


(56a) 


where 


- -uu 


A r — | 0 Ad e Bdte^pe 

_ B sa Cda B se C de A s 


B r a — 


Bdra 

0 

0 


B r 


0 

B d re 

0 


Cr = [0 0 C s 


(56b) 


where 6 a and 9 e are azimuth and elevation encoder readings, C pa and C pe are the output matrices for the 
azimuth and elevation pinion rates, and C a and C e are the output matrices for the azimuth and elevation 
encoders, respectively. 

Figure 17 shows the magnitude of the transfer function from the azimuth rate input r a to the azimuth 
encoder rate 9 a (solid line) and the magnitude of the transfer function from the elevation rate input r e to 
the elevation encoder rate 9 e (dashed line). The figure shows that the required identity relationship for 
low frequencies is not acquired. The magnitude of the transfer functions for frequencies less than 0.3 Hz 
is 0.74, below the required 1, due to inaccuracy in the model parameters (mainly in the hydraulic part). 
T his drawback can be removed by the experimental investigation of the parameters of the hydraulic 
drives, such as motors, valves, and lines. However, this inaccuracy is corrected by the position feedback 
loop, as will be shown later. The high-frequency peaks in azimuth and elevation (8 Hz in azimuth and 
20 Hz in elevation) are the gearbox resonances. 
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Fig. 16. Block diagram of the rate loop model. 



Fig. 17. Magnitude of the transfer function of the rate loop model. 


V. Position Loop Model 

The rate loop system with the proportional and integral (PI) controller is shown in Fig. 18, where e a 
and e e are the azimuth and elevation servo errors. For the series connection of the rate loop system and 
the controller, as in Fig. 18(a), define the state vector = [x ai x ei xj t ] with the new state variables 
x ei and x a i (integrals of the errors) such that 


ai — 1 



(57) 


The system output y is defined in Eq. (56a), the encoder output is 0 T = [6 a 0 e \, and the input is 

e T = [e a e e ]. The inputs to the rate loop systems are obtained from Fig. 18(a): 


X a — k pa &a T ^ 2 a X a i ) 


V e — k pe e e “I - ki e X e i 




(58) 


where k pe ,k ie ,k pa , and ki a are proportional and integral parameters of the controllers. Combining the 
equations for the rate loop system with Eqs. (57) and (58), one obtains 
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(a) 




Fig. 18. Position loop: (a) open and (b) closed. 


where 


Xq — A. 0 Xq “I - B 0 c 

9 = C 0 x 0 > 

y = Cx 0 


0 

0 

kia Bra 


0 

0 

kieB r e 


0 

0 

At 


Bo = 


I 

0 


0 

7 


kpaBra kp e I3re _ 


Co 


0 0 c a 
0 0 Ce 


C=[ 0 0 C r ] 


For the closed-loop system [see Fig. 18(b)], 


(59a) 


(59b) 


e = c — 6 (60) 

where c T = [c a c e ] is a command signal in azimuth, c a , and in elevation, c e . Introducing Eq. (60) to 
Eq. (59), one obtains 


x c i — x qi -f- B r c 


y = Cx c i 


(61a) 


where 
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A-ci — A 0 B 0 C c 


(61b) 


The simulations shown in Figs. 19 through 22 have two sets of assumptions: (1) a proportional gain 
of 1 in azimuth and elevation and an integral gain of 0.3 and (2) a proportional gain of 0.7 in azimuth 
and elevation and an integral gain of 0.2. The closed-loop transfer functions from azimuth command 
to azimuth encoder are shown in Fig. 19(a), and those from elevation command to elevation encoder 
are shown in Fig. 19(b). They show a bandwidth of 0.1 Hz. The cross-coupling transfer functions from 
azimuth command to elevation encoder and from elevation command to azimuth encoder are shown in 
Fig. 20. They show low-level cross-coupling. The closed-loop step responses from azimuth command 
to azimuth encoder are shown in Fig. 21(a), and those from elevation command to elevation encoder 
are shown in Fig. 21(b). They show a 20- to 30-percent overshoot and a 7- to 9-s settling time. The 
cross-coupling from azimuth step command to elevation encoder and from elevation step command to 
azimuth encoder is shown in Fig. 22. The cross-coupling is of the order 10 -3 . 
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Fig. 1 9. Magnitude of the transfer function of the 
closed-loop system: (a) azimuth command to 

azimuth encoder and (b) elevation command to 
elevation encoder. 
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Fig. 20. Magnitude of the transfer function of the 
closed-loop system: azimuth command to 

elevation encoder and elevation command to 
azimuth encoder. 



01 23456789 10 

TIME, s 


Fig. 21. Step responses of the closed-loop 
system: (a) azimuth command to azimuth 

encoder and (b) elevation command to elevation 
encoder. 





Fig. 22. Step responses of the closed-loop 
system: azimuth command to elevation encoder 
and elevation command to azimuth encoder. 


VI. Wind Disturbance Simulations 

Wind gust disturbances were modeled similarly to the DSS-13 antenna (see [13]) using the wind tunnel 
pressure distribution on the dish taken from Blaylock. 1 Their time history is generated using the wind 
Davenport spectrum (see [14] and [15]), determined for the Goldstone site. The simulations for the 
50 km/h wind gave the results listed in Table 4 and compared with the simulation results of the DSS-13 
antenna. The table shows that DSS 14 has better disturbance rejection properties (at the encoders) than 
has the DSS-13 antenna. 


Table 4. Servo errors in mdeg (3cr rms) 
for 50 km/h wind gusts. 


Drive 

Front wind 

Side wind 

Elevation, DSS 14 

2.6 

0.7 

Elevation, DSS 13 

14.6 

1.9 

Azimuth, DSS 14 

0.1 

2.1 

Azimuth, DSS 13 

0.5 

2.3 


VII. Conclusions 

An analytical model of the DSS-14 antenna has been developed. The rate loop model consists of the 
structural model (derived from the finite element model), gearbox model, hydraulic servo, and electronic 
boxes. The position loop was closed, and the time and frequency responses were simulated. The wind 
pointing errors of the DSS-14 antenna have been simulated. The model allows for detailed simulation of 
antenna dynamics and for modifications and improvements to the antenna control system. 

The simulations confirmed that the use of encoders located at drives limits the performance of the 
antenna (mainly by reducing its bandwidth to 0.1 Hz). The use of the master equatorial or new encoders 


1 R. B. Blaylock, “Aerodynamic Coefficients for Model of a Paraboloidal Reflector Directional Antenna Proposed for a 
JPL Advanced Antenna System,” JPL Interoffice Memorandum CP-6 (internal document), Jet Propulsion Laboratory, 
Pasadena, California, 1964. 
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located close to the axes of rotation of the antenna (similarly to the 34-m antennas) would allow expansion 
of the bandwidth to 0. 7-1.0 Hz. 


The antenna model needs further improvement. First, in this model, certain parameters of the hy- 
draulic drive are known with rather poor accuracy, and it influences the accuracy of the antenna model. 
It is essential to use experimental techniques to get more precise values of the parameters. Secondly, the 
RF pointing errors (in elevation and cross-elevation) of the antenna should be determined in order to 
evaluate the precision of the antenna pointing. 
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Appendix 

Transfer Function Derivation 


Each component of the electronics board is composed of operational amplifiers (opamps), resistors, 
and capacitors. The basic configuration of an inverting opamp circuit is shown in Fig. A-l. The “+” 
terminal of the opamp is grounded; thus, the “— ” terminal voltage is zero, called a virtual ground. In 
this situation, the currents ii and i 2 flowing through impedances Z\ and Z 2 are equal to 


} (A-l) 


and their sum is zero; that is, ii = — i 2 . Introducing them to Eq. (A-l) gives 






*2 


Vout 

z 2 


Vout — hv^n 


(A-2a) 


where 



(A-2b) 
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Fig. A-1 . Opamp circuit. 


I. Transfer Function Gf 0 

A schematic for the transfer function G to is shown in Fig. A-2(a), where the notation and the value 
of each element were taken from JPL Drawing 9479871D. 2 This schematic can be simplified to the one 
shown in Fig. A-2(b). In this figure, 

R s = {^R^q T T + R 59 ^) + Fi ’)2 — 49.5 kQ (A-3a) 

where R 5e = i ? 57 = R58 = R59 = 100 kSl, and Rq2 = 24.5 kfi; thus, R s = 49.5 kfi. The component Z\ is 

^1 = i?63 + = ^3 + i?64 - 91.1 m (A-3b) 

i + rt64^40 s 

where R$s = 40 kO, R $ 4 = 51.1 kf2, and C 40 = 0.15 fiF. The time constant R 64 C 40 — 0.0077 s is small, 
thus neglected. Denote Rmta — 9.7 kfl and R rn t e = 7.8 kfi the motor resistances in azimuth and elevation, 
respectively, and C 41 = 0.15 /jF. Then, the component Z 2 for the azimuth drive is as follows: 

Z 2 = 7 7 p rr,t V - - R mta = 9.7 kfi (A-4a) 


and for the elevation drive, 


Z 2 = i T ~ ^ = 7 ‘ 8 (A-4b) 

1 -r -fL 7 n f :e 04 xS 

The time constants R m taC 4 \ = 0.0015 s and R m teC 41 = 0.0012 s are of the order 10~ 3 s, thus considered 
small, and neglected. 

The transfer function G t0 for azimuth is 


Gto — 


R P 

R s + Rp 


8800 

49, 500 + 8800 


0.151 


(A-5a) 


and for elevation, it is 


2 JPL Drawing 9479871D (internal document), Jet Propulsion Laboratory, Pasadena, California. 
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(A-5b) 



Fig. A-2. Schematic for the transfer function G/q: (a) full and (b) simplified. 


II. Transfer Functions G r i and G r 2 

The transfer functions G r \ and G r 2 are determined simultaneously. Their schematic is given in 
Fig. A-3(a), the parameters parameters of which are 

R i5 = 750 kft 

R 50 — 100 kfi 

i?5i = 12.1 kfl 

R 52 = 442 kfl 

R 53 = 442 kO 

Rq5 — 909 kO 

i?66 = 90.9 kfl 

C 31 — 1 fJ-F 

C 42 = 0.1 fiF 

The schematic from Fig. A-3(a) can be transformed to the form shown in Fig. A-3(b). The value of Z 3 
is as follows: 

Zz = i?66 + 1 p 65 /^ — -^66 + #65 = 10 6 kfl (A-7) 

1 + i<65W2S V 
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In this variable, the small time constant ResC^ = 0.0909 s was ignored. 


The value Z 4 is obtained as 


Z 4 = R 53 4- 


Rp 

1 + R 0 C 31 s 


- 4.65 x 10 6 


1 + 0.400s 
1 + 4.205s 


(A-8) 


where R 0 = (R50 + R52 + R50R52)/ R51 = 4205 kO. 

Having determined Z 3 and Z 4 , the transfer functions G r 1 (from v r to v s ) and G r 2 (from v to to v s ) are 
obtained: 


G r i=~ = 6.20 G a 

“15 

G T 2 = — ~ = —4.65 G 0 
63 


> 


(A- 9 ) 


where 


1 + 0.400s 
0 ~ 1 + 4.205s 


(A-10) 


is the transfer function of a lag compensator. 




Fig. A-3. Schematic for the transfer functions and G^: (a) full and (b) simplified. 


III. Transfer Function G s 

The transfer function G a is determined from the schematic in Fig. A-4(a), and is shown in compact 
form in Fig. A-4(b). For this schematic, Ri 3 = 100 kO, R 36 = 10 kQ, R43 = 24.9 kfl, and Gis = 0.1 /iF; 
therefore, one obtains 


^5 


R36 

1 + R36G18S 


“ R 36 = 10 kO 


(A-ll) 


where i2 36 Ci8 = 0.0015 s = 0. Since Vi — v s Z 5 /R 43 , and i s = i>i (Z 5 1 + R 1 ), thus, 
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